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Abstract 

We analyze and interpret recent optical experiments with semiconductor quantum 
dots. We derive a quantitative relation between the amount of information trans- 
ferred into the environment and the optical polarization that may be observed in a 
spectroscopy experiment. 
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1 Introduction 



The fundamental difficulty in sustaining the coherence of a quantum system is 
its interaction with the surrounding world. In the course of joint evolution of 
the system and its environment, this interaction establishes phase correlations 
between the former and a macroscopic number of degrees of freedom of the lat- 
ter. For most states, such correlations perturb and eventually erase the internal 
phase information of the quantum system turning a non-classical superposi- 
tion state into a mixture of a small subset of classicaly allowed states. This 
effect is known as decoherence or dephasing and seems to be one of the most 
fundamental aspects of the quantum theory [1,2]. It is also of importance for 
certain practical tasks, like quantum computing [3], where maintaining system 
coherence over many control operations is indispensable. 

Once the quantum correlations (entanglement) between the system and its 
environment have been created, a measurement on the environment may, in 
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principle, yield information on the system. Thus, one may say that a cer- 
tain amount of information, allowing one to determine the system state, has 
been transferred to the environment. Since both information transfer and de- 
coherence result from the same underlying physical process of building up 
correlations, a certain relation between the amount of information transfer 
and the phase coherence retained in the system should be expected. The most 
celebrated textbook example of such a relation is the vanishing of interference 
fringes in a double-slit experiment with single particles whenever one tries to 
establish through which slit the particle passes. Because of this historical rela- 
tion to the interference of two spatial paths, the knowledge of the system state 
is customarily referred to as "which path" (or welcher Weg) information. 

The decoherence effect due to the transfer of information to the environment 
appears always when a quantum system is coupled to the surrounding world. 
In this paper we revisit one specific example of this effect: the decay of coherent 
optical polarization in a quantum dot due to the dephasing of confined carrier 
states, which we will interprete as a result of information leakage from the 
carrier subsystem to its environment (phonon modes). After introducing the 
system and defining its model (Sec. 2) we will derive the evolution of the 
confined carriers coupled to phonons using an algebraic technique based on 
Weyl operators (Sec. 3). Next, we will discuss a quantitative measure of the 
distinguishability of quantum states due to the information contained in the 
environment (Sec. 4). Finally, in Sec. 5, we will use the facts introduced in the 
preceding sections to show that the pure dephasing effect observed as the decay 
of the coherent response from an optically excited semiconductor quantum 
dot (or an ensemble thereof) [4] may be interpreted as a result of which path 
dephasing. We derive a quantitative complementarity relation between the 
observed coherent optical response and the amount of information transferred 
to the environment. 



2 The system and the model 

In this paper, we consider the simplest version of a time-resolved optical ex- 
periment performed on a single quantum dot (for practical reasons, actual 
experiments are often performed on ensambles of QDs using nonlinear tech- 
niques [4]). A very short laser pulse prepares the system in a certain super- 
position (dependent on the pulse phase and intensity) of the ground state (no 
exciton, denoted |0)) and the single-exciton state (denoted |1)). By very short 
we mean a pulse much shorter than the time scales of phonon dynamics, so 
that the preparation of the initial state may be considered instantaneous. This 
corresponds to the actual experimental situation with pulse durations of order 
of 100 fs [4]. On the other hand, the pulse is long enough to assure a relatively 
narrow spectrum and prevent the population of higher confined levels. Let us 



2 



restrict (for simplicity) to the equal superposition state 



« = (i) 

In such a state, the inter-band component of the electric dipole moment has a 
non- vanishing average value oscillating at an optical frequency (hence referred 
to as optical polarization) which leads to the emission of coherent electromag- 
netic radiation with an amplitude proportional to the oscillating dipole mo- 
ment. In an unperturbed system (e.g., in an atom), the radiation would be 
emitted over time of the order of the lifetime of the superposition state, i.e., 
until the system relaxes to the ground state due to radiative energy loss. 

In a semiconductor structure an additional effect, related to carrier-phonon 
coupling, appears on a time scale much shorter than the lifetime of the state. 
Due to the interactions between confined carriers and lattice ions, the ground 
state of the lattice in the presence of a charge distribution is different than 
in its absence. As a result, after the creation of a confined exciton the lattice 
relaxes to a new equilibrium, which is accompanied by the emission of phonon 
wave packets [5,6] that form a trace in the macroscopic crystal distinguishing 
the exciton state from an empty dot. As we will discuss below, this information 
broadcast via emitted phonons leads to a decay of the coherence of the super- 
position state although the average occupations of the system states remain 
unaffected (hence the process is referred to as pure dephasing) . Since coherent 
dipole radiation requires well-defined phase relations between the components 
of a quantum superposition, the amplitude of this radiation, measured in the 
experiment, gives access to the coherence properties of the quantum state of 
confined carriers itself. The dephasing of the qantum superposition is there- 
fore directly translated into the decay of coherent optical radiation from the 
system. 

For the carrier-phonon system in a semiconductor quantum dot, a microscopic 
model of the dephasing effect exists (involving the interaction with acoustic 
phonons), which reproduces experimental data very well [7,8] and may serve 
as a reliable starting point to describe the evolution of the combined system 
of confined carriers and lattice modes. 

The Hamiltonian of the system is 

H = e|l)(l| + H ph + |1)(1| J2(f*kh + f h bl), (2) 

fc 

where the first term describes the energy of the confined exciton (e is the energy 
difference between the states without phonon corrections), H ph = Sfc^fc&fc&fc 
is the Hamiltonian of the phonon subsystem and the third term describes the 
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interaction. Carrier-phonon interaction constants in (2) are given by 



/ n h f°° 

fk={aC ~ ah) ^^cJ- o d3rr{r)e ~ lk ' T ' P{r) ' (3) 

and describe the deformation potential coupling between the carriers and the 
lattice modes, which is the dominating mechanism under the assumed optical 
driving conditions [9]. Here g is the crystal density, Vn is the normalization 
volume of the phonon system, = ck is the frequency of the phonon mode 
with the wave vector k (c is the speed of longitudinal sound), and b\,, bj~ 
are phonon creation and annihilation operators. The exciton wave function 
is modelled as a product of two identical single-particle wave functions ip(r e ) 
and ^>(rh), corresponding to the electron and hole, respectively. 

In our calculations we use typical parameters for a self-assembled InAs/GaAs 
structure: single particle wave functions tp(r) modelled by Gaussians with 
4 nm width in the xy plane and 1 nm along z, the deformation potential 
difference a e — a h = 9.5 eV, crystal density g = 5300 kg/m 3 , and the speed of 
longitudinal sound c = 5150 m/s. 



3 Phonon-induced pure dephasing of optical polarization 



In this Section we describe the dephasing of a confined exciton state within the 
exactly solvable model of interaction with the environment presented above. 
After describing the time evolution of the system we find the degree of co- 
herence remaining in the system as manifested by the amplitude of the ex- 
perimentally measurable coherent dipole radiation. In this way we reproduce 
the recent theoretical description [9] using a simple algebraic method which 
provides the complete density matrix of the carrier subsystem, necessary for 
the discussion to be presented in Sec. 5. 

The carrier-phonon interaction term in Eq. (2) is linear in phonon operators 
and describes a shift of the lattice equilibrium induced by the presence of a 
charge distribution in the dot. The stationary state of the system corresponds 
to the exciton and the surrounding coherent cloud of phonons representing 
the lattice distortion to the new equilibrium. The transformation that creates 
the coherent cloud is the shift wb^w^ = bj~ — /*./ (fkVk), generated by the Weyl 
operator (see Appendix) 



w = exp 



(4) 



A straightforward calculation shows that the Hamiltonian (2) is diagonalized 
by the unitary transformation W = |0)(0| <g> I + |1)(1| <S> w, where I is the 
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identity operator and the tensor product refers to the carrier subsystem (first 
component) and its phonon environmnent (second component). As a result 
one gets 

H = WHW ] = E\l)(l\ + H ph , 
where E = e - £ fc \f k \ 2 /(hw k ). 

We assume that at the beginning (t = 0) the state of the whole system is 
Po = (IV'oXV'oI) ® Pe, where pe is the density matrix of the phonon subsystem 
(environment) at thermal equilibrium and \ip ) is the equal superposition state 
(1) prepared by a properly chosen ultrashort pulse. 

The evolution operator U (t) = e~ lHt / h may be written as 

U(t) = W^WU(t)W^W = W ] U{t)W = W ] U{t)WU\t)U{t) = W^W(t)U(t), 

where U(t) = e~ iSt / h and W(t) = U(t)WW(t). Since U(t) is diagonal the 
explicit form of W(t) may easily be found and one gets 



U(t) = [|0><0 1 ® I + |1)(1| ® w*w(t)] U(t), 
where w(t) = e~ lH P bt ^ h we lHpbt ^ h . 



(5) 



Using the evolution operator in the form (5) the system state at a time t may 
be written as 



1 



Pe 



z iEt/n p E wt(t)w 



p(t) = - . . . 

2 \ e - lEt l h w^w{t)pE W^w(t)pEVj\t)w J 



(6) 



where we used the tensor product notation in which an operator A is expanded 
as A = J2mn \ m )( n \ ® A mn with a set of operators A mn acting on the second 
subsystem, and written in the matrix form with respect to the first subsystem. 
The density matrix for the carrier subsystem is obtained by tracing out the 
phonon degrees of freedom, i.e., ps = Trs p. Hence, 



Ps(t) 



1 e iEt l h {w\t)w) ^ 

e- iEt ' h (w^w(t)) 1 ) 



(7) 



where the average may be calculated by first using the multiplication rule for 
Weyl operators [Eq. (A. 3)] to combine and w(t) and then Eq. (A. 5). The 
result is 



[w 





fk 


exp |-2 





[i sin uo k t + (1 - cos u k t) (2n k + !)]>, 



where n k are bosonic equilibrium occupation numbers. 
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Fig. 1. Decay of the coherent radiation from a confined exciton at various temper- 
atures, as shown. 

The emitted coherent dipole radiation is proportional to the non-diagonal 
element of the density matrix p s (t) and its amplitude is 



P(t) = P \{w\t)w)\. 



(8) 



In Fig. 1 we show the normalized polarization amplitude P(t)/P (first derived 
in Ref. [9]). The interaction with the macroscopic crystal environment leads to 
a reduction of coherent radiation due to pure dephasing of the exciton state, 
reflected by the reduced value of the non-diagonal element of the density 
matrix pg. At t — one has (w^(t)w) = 1, while at large values of t, cosujkt 
oscillates very quickly as a function of k and averages to (see also Ref. 
[9]). Thus, for long times, the polarization amplitude tends to a temperature- 
dependent finite value 



P(t) ^P exp 



fc 



fk 



(2n fc + 1) 



<Po- 



This partial decay of coherence is a characteristic feature of short-time de- 
phasing for carrier-phonon couplings encountered in real systems [4]. 

On the other hand, it is clear from the diagonal elements in Eq. (6) that the two 
carrier states |0) and |1) are accompanied by different states of the phonon en- 
vironment so that correlations between the two subsystems are present. These 
correlations may be related to the phonon wave packets that are emitted into 
the bulk of the crystal after ultrafast optical excitation [6,5]. These wave pack- 
ets traveling away from the QD region carry the information about the system 
state into the environment. We will see that dephasing may be quantitatively 
related to the amount of these correlations, i.e., to the information on the 
exciton state extracted by its phonon environment. First, however, we need a 
quantitative measure for this information. 



4 Distinguishability of quantum states 



For completeness, in this Section we review the definition of a measure for 
the quantity of which path information on a quantum state of the system 
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contained in its environment [10]. 

Qualitatively, the definition is as follows. Let us consider a two-level system S 
with the Hilbert space spanned by the states |0), |1), interacting with its envi- 
ronment E. A which path measurement of the system S, e.g., a measurement 
of the observable |0)(0|, determines the system state. Given the state of the 
total system S+E, how well can we predict the outcome of this measurement 
based on a previous measurement on the environment? The measure of distin- 
guishability is based on the probability p of a correct guess for the best possible 
choice of the measurement on the environment. If the systems are completely 
uncorrelated then any measurement on the environment is of no help and the 
probability of correct prediction remains equal to 1/2. If the systems are in 
a maximally entangled state like |0)|£b) + with (E \Ei) = then the 

result of the measurement of the environment observable \E )(E \ completely 
determines the state of the system S, so that p — 1. The distinguishability 
measure is defined as 

V = 2 (p-\) ( 9 ) 
and changes correspondingly from to 1. 

In general, for a given state of the total system p and a result e of the en- 
vironment measurement, the optimal prediction procedure [11] is the usual 
statistical inference rule: one calculates the probabilities 

p(i,e) = Tr[(\i)(i\®\e)(e\)p\ (10) 

of finding the system S in the state \i) simultaneously with the environment 
in the state |e) associated with the measurement outcome e. Then, given the 
outcome e, one opts for the state % for which this probability is larger. The 
prediction is correct if the total system is indeed in this more probable state, 
whichever it is. Allowing for all possible outcomes e, the total likelihood of 
guessing right is equal to the probability that the system is in any of the "more 
probable states" , i.e, 

p = ]T max [p(0, e),p(l, e)] = W [p(°> e ) + p( x > e )] + \ E b(0, e ) " p(!> e ) I ■ 

e ^ e ^ e 

The first term contains the sum of probabilities for any state and is therefore 
equal to 1/2. Writing the full density matrix (for equal probabilities of the two 
system states) in the form p = (1/2) ® Piji where poo and pu (but 

neither p i nor Pio) are density matrices, and using Eq. (10) we can write 

P = \ + \H I < e l/°oo -Pn|e)| • 

This formula gives the probability of correct prediction for a fixed measure- 
ment on the environment, i.e. for a specific basis of environment eigenstates 
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(|e)}. In order to measure the amount of information inherently contained in 
the environment, independently of the possibly poor choice of the measure- 
ment basis, one maximizes this quantity with respect to all possible observ- 
ables that can be measured, i.e. to all possible complete eigensystems {|e)}. 
To find the optimal (upper bound) value for p, we split the hermitian opera- 
tor Ap = poo _ Pn into a positive and negative part. Let Api and \l) be the 
eigenvalues and the corresponding eigenstates of Ap. Let L + , L_ be the sets of 
quantum numbers / for which Api is nonnegative and negative, respectively. 
Then Ap = Ap^— Ap<~), with the positive operators ApW = J2ieL± ^Pi\ 
Since these operators are positive, one has 



(e|ApW|e) - (e|Ap^|e) < (e\Ap (+ >\e) + (e\Ap^>\e) = (e ||Ap|| e 
where | Ap| is the modulus of the operator Ap, defined as 
|Ap|=^|Ap ; ||/)(/|=ApW + Ap(-). 



Hence, 

\^\(e\Ap\e)\<\Tr\Ap\ = ±D(p 00:Pll ) : 

where D(p, p') is known as the trace distance [3] between the density matrices 
p and p' . On the other hand, the equality is clearly attained if the eigensystem 
(|e)} coincides with the eigenstates of Ap, i.e., when the observable measured 
on the environment is Ap itself. Thus, using the definition (9), we finally get 

V = D( Poo , Pll ). (11) 



5 The complementarity relation 



In this Section we derive a quantitative complementarity relation between 
the degree of coherence in the quantum dot system, as manifesteed by the 
amplitude of coherent radiation described in Sec. 3, and the amount of which 
path information transferred to the environment, as defined in Sec. 4. This 
relation is analogous to the visibility-distinguishability relation in the double- 
slit setup [10,12]. 

For the carrier-phonon state of Eq. (6) the distinguishability of carrier states 
due to the correlations with the phonon environment is, from Eq. (11), 

V{t) = D(p E ,w ] w(t)pW(t)w). 
One has, in general [3], D 2 (p, p') < 1-F 2 (p, p'), where F(p, p') = Tr yj p 1 / 2 p 1 p 1 / 2 



8 



is the fidelity measure of the distance between the states p and p' . Hence 

2 



V\t) < 1 - 



Tr \J p 1 J ] 2 w^w(t)p E w^(t)wp 1 l 



/2 



= 1 



Tr 



1/2 + (4 x i/ 



1/2 



(12) 



Since Tr ] A\ > | Tr A\ we may write 

Tr Ip^V^OAe 2 ! > I Tr(j> E w*w(t))\ = \(w^w(t))\ = P(t)/P , (13) 
where we used Eq. (8). Combining Eqs. (12) and (13) leads to the relation 

P(f)" 12 



+ V\t) < 1, 



(14) 



which shows that the relative decay of polarization is related to the which 
path information transfer to the environment. Note that for pure states pe, 
equality holds in Eqs. (12) and (13), so that the relation (14) also turns into 
an equality. 



6 Conclusion 



In this paper, we have interpreted the dephasing effects observed in recent 
optical experiments on semiconductor quantum dots in terms of "which path" 
information transfer to the environment. Coherence properties of carriers con- 
fined in a QD may be tested by detecting coherent radiation emitted from the 
system as a result of optical polarization induced by optical excitation. Due to 
the coupling between the system and its phonon environment, in the course of 
quantum evolution lattice modes get excited, which leads to the dephasing of 
the carrier state manifested by a decrease of the amplitude of this radiation. 
We have derived an inequality between the remaining amplitude of emitted 
radiation and the amount of information carried to the bulk of the macro- 
scopic crystal by lattice excitations. From this point of view the dephasing of 
optical polarization is of the same nature as the well-known disappearance of 
interference fringes in a double-slit (Young) experiment due to the which path 
knowledge. 

The complementarity relation (14) is analogous to similar relations between 
the distinguishability of paths and the visibility of interference fringes in the 
double-slit setup [10,12]. In fact, an interference experiment analogous to the 
double-slit experiment has been performed using optically excited states of 
semiconductor quantum dots [13], where the two states (analogous to paths) 
are the absence or presence of a single, optically created exciton, confined by 
the binding potential of the QD. It may be shown that the visibility of inter- 
ference fringes in these time-domain interference experiments is governed by 
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the same non-diagonal element of the reduced density matrix as the coherent 
polarization in our discussion. 



Our discussion confirms that incompatibility of the quantum behavior and 
which path information is a general feature of the quantum world. We have 
shown that a manifestation of this general property may be identified in an 
experimental situation that is extremely different from the historical concept 
of two spatial paths in a double-slit experiment. 
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Appendix 

A Rules for Weyl operators 

Applying the Weyl operator (4) to the phonon annihilation operator b k and 
expanding the result we get [14] 



where S = J2k(gkb k -g k b k ), g k = f k /(huj k ). Since [S,b k } = -g k and [S, [S, b k ]]=0 
(as all the higher commutators), we obtain the shift effect of the Weyl operator, 
wb k w^ = b k — g k . 

Multiplying Weyl operators is done using the Baker-Hausdorff formula for 
separating operator exponents [14] 



where A and B are some operators and [A, [A, B]] = [B, [A, B]] = 0. Some 
simple algebra using Eq. (A.l) shows that the product of two Weyl operators, 





(A.2) 



i — 1, 2, is equal to 



1 ( (1)* (2) (1) (2)*\ 

Wl w 2 = w 3 exp -- }_^[g k g k - g k 'g k ) 



k 



(A.3) 
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where w 3 is given by (A. 2) with g k 3) = g k l) + g k 2) . 

Now it is possible to calculate the equilibrium average of a Weyl operator 
(following Ref. [14]), (w) = Tt(pew), where ps is the density matrix of the 
phonon reservoir at thermal equilibrium. Since pe = e~^ u ^ j Tr e _/3Hph , (3 = 
l/(k B T), and H ph = Y.k^kb k b k , we may write 

Tr (e~^fe^ fe6 fe 6fe e£fe (9fe6 fe~ 9 fe 6fe) l 
( w ) = ^ L 

and separate the thermal equilibrium into 

(w) = U(w) k , (A.4) 
k 

where 

{w)k = f ( A - 5 ) 

Using Eq. (A.l) and keeping in mind that Tr(e~ l3hjJhb i bh ) = (1 - e-^y 1 , we 
can transform Eq. (A. 5) into 

( w ) h = e ^\l - e-P**") Tr(e-^ fe6 fe 6fe e 9fe6 fee-^ 6fe ). (A.6) 



For the trace one has 

t t 00 t 

Tr^-Z^^e^e-s^) = £ e -^ fcm (m|e fffc6 *e-^ bfc |m), (A.7) 

m=0 

where |m) denotes a state with m excitations. The exponents are expanded in 
a power series, 

oo (_„*\l 

e -g* k b k | m j = ^ L^L &fc | m ). (A.8) 
i=o 

Since 6' fc |m) = [m!/(m — Z)!] 1,/2 |m — /) for / < m and b l k \m) = for Z > m, 
Eq. (A.8) transforms into 



-a* k K 



m ( „*\l 



e 7TZ ) 



E 



7! 

2=0 l ' 



to! 



1/2 



(to-/)! 

Treating the other exponent in the same manner we end up with 



m-l). (A.9) 



. m (_\„ \2\l m \ 



11 



that is, a Laguerre polynomial of order m, L m (\gj~\ 2 ). Using the fact that 

oo 
m=0 

with a specific z = e~ l3hulk [so that = z/(l — z)], and inserting the result 
into Eq. (A. 6), and then into Eq. (A. 4), we get the thermal average of the 
Weyl operator in the form 

( W ) = e -5Ej9 fe | 2 (2« fe +l). (A.ll) 
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